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Abstract 

We study the dynamical response of the 16 O nucleus to an incident antiproton using the Giessen 
Boltzmann-Uehling-Uhlenbeck microscopic transport model with relativistic mean fields. A special 
emphasis is put on the possibility of a dynamical compression of the nucleus induced by the 
moving antiproton. Realistic antibaryon coupling constants to the mean meson fields are chosen 
in accordance with empirical data. Our calculations show that an antiproton embedded in the 
nuclear interior with momentum less than the nucleon Fermi momentum may create a locally 
compressed zone in the nucleus with a maximum density of about twice the nuclear saturation 
density. To evaluate the probability of the nuclear compression in high-energy p-nucleus collisions, 
we adopt a two-stage scheme. This scheme takes into account the antiproton deceleration due to 
the cascade of pN rescatterings inside the nucleus (first stage) as well as the nuclear compression by 
the slow antiproton before its annihilation (second stage). With our standard model parameters, 
the fraction of p annihilation events in the compressed zone is about 10~ 5 for p 16 collisions at 
Plab = 3 — 10 GeV/c. Finally, possible experimental triggers aimed at selecting such events are 
discussed. 
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I. INTRODUCTION 



The production of compressed nuclear matter in laboratory is one of the most important 
achievements of heavy-ion physics during last decades. Heavy-ion collision experiments open 
the possibility to study new phases of matter, such as e.g. a quark-gluon plasma [l, 2] (see 
also [3] for a recent review). In a heavy-ion collision, compression is accompanied by the 
strong heating of matter by the shock wave mechanism [4j]. However, very little is known 
about possible compressional effects induced by a slowly moving or even stopped hadron in 
a nucleus. In this case, the compression is associated with the enhanced concentration of 
nucleons around the hadron, provided its interaction with nucleons is sufficiently attractive. 
Several examples of such systems are under discussion, but their existence is still an open 
question. 

The most famous example of strongly-bound hadron-nucleus systems are A-hypernuclei. 
By measuring the E2(5/2 + — > l/2 + ) transition in the 7 A Li hypernucleus, the shrinkage of 
the 6 Li core size by ~ 19% has been found experimentally in Ref. |5|. The hypothetical 
multi-strange nuclei composed of several hyperons (A, 5°, S~) and nucleons (see 6|, |7|] and 
references therein) could be selfbound and have an enhanced baryon density. 

As proposed in Ref. 8j on the basis of a phenomenological potential model, hypothetical 
A'-nuclei could be long-lived and strongly bound compact systems with the nucleon density 
reaching almost 10po, where po is the normal nuclear matter density. However, as stated 
in the recent work by Hayano et al [9j, the measurement of 2p shift in the 3d-2p transition 
in kaonic 4 He, eliminating a long-standing discrepancy between the standard theory and 
experiment, poses severe limitations on superstrong potentials with such high densities. 
Moderate compressional (core polarization) effects, up to about 2p n , have been found within 
the relativistic mean field calculations of A'-nuclei done in Ref. [lOj. However, the depths 



of 100-200 MeV for the real part of t 



re A'-nucleus potential at po, which follow from the 



phenomenological models of Refs. (si, 10], are disfavoured by the chiral SU(3) models [llL 
131 ] . Nevertheless, the existence of strongly-bound A^-nuclear systems largely related to the 
nature of A(1405) is still under theoretical debates Q, Q and experimental studies 

Recently strong compressional effects have been predicted in the strongly bound p-nucleus 
systems Q] in the case of a deep real part of an antiproton optical potential, Re(V r opt ) < 
— (150 — 200) MeV at po- We remark, however, that the antibaryon optical potentials in the 



nuclear interior are largely unknown and their study requires more efforts 19M25 ] . 



In the present work, we extend our previous study of the dynamical compression induced 



by a stopped antiproton 



26] to the case of a moving p. The calculations are based on 



the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model 271 ] . First, we study 
kinematical and geometrical conditions at which an antiproton can generate the increase 
of nucleon density. Second, by performing the transport simulations of p-nucleus collisions 
we evaluate the actual probability of p-annihilation in the compressed zone for the beam 
momenta of 0.3-10 GeV/c, relevant for future antiproton beams at the Facility for Antiproton 
and Ion Research (FAIR) in Darmstadt. Finally, we study possible triggering schemes which 
can be used to select the events with p-annihilation in the compressed nuclear environment. 



We 



13, 



lave chosen the O nucleus as a target. This is motivated by our earlier observation 



18 



261 ] that the compressional effects associated with p are more pronounced in light 



nuclei. 

The paper is organized as follows: In Sect. [Til the description of a calculational procedure 
is given. Then, in Sect. 1111} we study the dynamical patterns of the nuclear compression by an 
antiproton initialized at different momenta and positions inside a nucleus. Sect. IIVI contains 
our results on the probabilities of a p- annihilation in the compressed zone for energetic p- 
nucleus collisions. In Sect. HV] we also discuss possible triggers based on the fast proton 
emission and on the measurement of the energy deposition. We analyse the influence of the 
possible in-medium modifications of the p-annihilation rate and of the different antiproton 
mean field parameters on our results in Sect. |V] Summary and outlook are given in Sect. IVII 



II. THE CALCULATIONAL PROCEDURE 



In our calculations, we apply the GiBUU model 27j. This model solves the coupled set of 
semiclassical kinetic equations for various hadronic species: nucleons, antinucleons, mesons, 
baryonic resonances and their corresponding antiparticles. We use the relativistic mean field 



mode of calculations 



22. 



26 



28 



29] which provides a simple and natural description of both 



baryonic and antibaryonic mean fields by using the same Lagrangian. The kinetic equation 
for the hadron of the type j (j = p, n, A ++ , A + , A , A - ,..., 7r + , 7r°, ir~,... and respective 



antibaryons) reads as 

R / . F)m*\ f) 1 

/ i (x,p*) = J i [{/}] , (1) 



1 

Po 



P 1" + 1 



where fj(x,p*) is the phase-space density of the j-th type particles, p* is the kinetic four- 
momentum (j»q = ^p* 2 + (m;-) 2 ), m* is the effective mass, and is the field tensor. The 
l.h.s. of Eq. (TI]) describes the evolution of the phase-space density fj(x,p*) under the 
influence of the mean mesonic fields. The r.h.s. of Eq. (TjQ) is the collision integral 
describing the change of the phase-space density due to the particle-particle collisions and 
resonance decays. 

The kinetic equations <Q are solved by applying the standard test particle technique in 
the parallel ensemble mode. The phase-space densities are represented by the set of the 
point-like test particles: 

/,(^P*) = 7T7 L E 8{T ~ T t (t)) 5( P * - P *(0) , (2) 
y?' jv ens j = i 

where Nj is the number of physical particles of the type j, N ens is the number of parallel 
ensembles, and gj is the spin degeneracy. The test-particle representation (J2]) provides a 
simple solution of the kinetic equations (JT|) in terms of the Hamiltonian-like equations for 
the centroids (i-j(t), p*(t)) (c.f. Eqs.(2),(3) in Ref. |22|). The collision integral is simulated 
with the help of a usual geometrical collision criterion (c.f. Ref. 28]). 

The mean mesonic fields are determined from the nonlinear Klein-Gordon-like equations 
with the source terms given by the particle densities and currents. Therefore, in order to 
provide a smooth coordinate dependence of the mean mesonic fields, the coordinate space 
^-functions in the r.h.s. of Eq. ([2]) are replaced by the Gaussians of the width L ~ 0.5 — 1 
fm in actual calculations. Then, e.g. the coordinate space density and the scalar density of 
the j-th type hadrons are computed as 

« w - I *m*, p*) - "f exp . (3) 

<*<*> - A I - wWf iM- 1 ^} ' (4) 

We are interested, in particular, in the values of the nucleon density p = p p + p n . The 
antiproton density p p and the nucleon scalar density p$ = Ps P + Psn are also used in the 
present analysis. 



The width L in Eqs.(j3]),(j4]) is a pure numerical parameter of the GiBUU model. Its value 



is correlated with 
step size (c.f. 



26 



;he number of parallel ensembles and is set equal to the coordinate grid 



28|). The physical results do not depend on L, provided that it is small 



enough to resolve the physical nonuniformities of the system. In the present calculations we 
use the value L = 0.5 fm from our earlier work [26j, where we have also studied the influence 
of L on the compression dynamics. 

For the nucleon mean field we apply the nonlinear Walecka model. The nucleon-meson 
coupling constants and the parameters of the cr-field self-interactions are taken from the NL3 
parameterization [3^]. This parameterization provides the nuclear matter incompressibility 
K = 272 MeV and the nucleon effective mass m* N = 0.6rriN at po — 0.148 fm -3 . Within the 
NL3 set of parameters, the binding energies, charge and neutron radii of spherical nuclei as 
well as deformation properties of some rare-earth and actinide nuclei have been described 
quite well 30] . The isoscalar monopole resonance energies in heavy spherical nuclei are also 
reproduced by this set of parameters [30]. 

The antinucleon-meson coupling constants are more uncertain. As it is well known, the 
G-parity transformation of Walecka-type Lagrangians results in too deep antiproton optical 



potentials. Therefore, following Refs. 



18 



19 



22 



26], we introduce a common reduction 



factor £ < 1 for the antinucleon coupling constants to the a-, u- and p- mesons as given 
by the G-parity transformation. Below, if it is not explicitly stated otherwise, we use the 
value £ = 0.22 obtained in [22] from the best fit of p-absorption cross sections on nuclei at 
the beam momenta below 1 GeV/c. The corresponding real part of an antiproton optical 
potential is about -150 MeV in the nuclear centre, which is somewhat deeper than the real 
part derived from the most recent p-atomic calculations [3], however, within the commonly 
accepted uncertainty interval^. 

Due to a big annihilation cross section, in majority of events, an antiproton colliding with 
a nucleus will annihilate already on peripheral nucleons. However, as argued in Ref. fisl ]. 
compressional effects are expected only in events when the antiproton penetrates deep to 
the nuclear interior and stops there due to (in)elastic collisions with nucleons. Such events 
are presumably quite rare and their study requires to go beyond the ensemble-averaged 
description provided by the kinetic mean field theory. The Quantum Molecular Dynamics 



3JJ or Antisymmetrized Mo 



1 For detailed discussion, see 



ecular Dynamics [32| models seem to be better theoretical tools 



22f and references therein 
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for studying such rare events. However, to our knowledge, at present there exists no version 
of a molecular dynamics model which incorporates all relevant antibaryon-baryon collision 
channels and relativistic potentials. 

In the present work, we treat compressional effects in a p-nucleus collision perturbatively. 
It means, that the influence of the compressional response of a nucleus on the deceleration 
process and eventual annihilation of an antiproton is neglected. Thus, the collisional dy- 
namics of the incident antiproton is simulated within standard GiBUU until its annihilation. 
We assume further, that the position and momentum of the p at the beginning of compres- 
sion process are not much different from those at its annihilation point. Then we study 
the compressional response of the nucleus to slow antiprotons and evaluate their survival 
probability. 

Therefore, we adopt a two-stage calculational scheme: On the first stage, an antiproton 
penetrates into the nucleus while experiencing one or more rescatterings on nucleons. We 
describe this process by the standard GiBUU simulation in the parallel ensemble mode with 
iVens = 1000 parallel ensembles. Each parallel ensemble is considered as one event. For each 
impact parameter, N IWS = 100 simulation runs have been done which gives iV ens A^ uns = 10 5 
events per impact parameter. We have chosen 32 impact parameters b = 0.25, 0.50, 8 fm 
for the p 16 system. Final results are impact-parameter weighted. Since the incoming p can 
be transformed to another antibaryons, e.g. n or A, we consider below the antibaryon anni- 
hilation in general. The coordinates and the kinetic three-momenta of an antibaryon 
just before the annihilation or, for events without annihilation, at the end of the computa- 
tional time (40 fm/c) have been determined and stored for every event. In the following, we 
always deal with the kinetic three-momenta of particles, but omit the word "kinetic" and 
the star symbol for brevitj^. Due to the averaging of the mean field over parallel ensembles, 
the compressional effects are practically unnoticeable in the standard GiBUU calculation, 
because rare events with a deep penetration of p into the nucleus are diluted with the ma- 
jority of events when the antiproton annihilates on the nuclear periphery. This is why we 
use the coordinates and momenta of the antibaryon obtained on the first stage as an input 



for another simulation based on the GiBUU model 



26 



331 ] . Thus, on the second stage, an 



antiproton is initialized inside the nucleus at the phase-space point (rg, p#) using the Gaus- 
2 In fact, if the collective motion of nuclear matter is negligible, e.g. when a fast hadron passes through 

the undisturbed nuclear target, the space components of the canonical and kinetic four-momenta are 

practically the same. 



sian distribution in coordinate space and the sharp-peaked distribution in momentum space. 
By doing so we neglect the possibility that the annihilating antibaryon can be different from 
the antiproton. This is, however, not important in view that the mean field contributions, 
apart from small isospin and Coulomb effects, are the same for all antibaryons in our model. 
The corresponding phase-space density of an antiproton is written as (H = c = 1): 

/p " (r ' P) = (2^7V3 6XP 1 2aTj ^- 5(P - P * } ' (5) 

where gp = 2 is the spin degeneracy of an antiproton and a r is the width of the coordinate 
space Gaussian. Please, notice, that the quantity o> is a physical parameter of our model, 
while the quantity L in Eqs. (J3|) , (HJ) is pure technical and should not be misidentified with a r . 
We stress that now the antiproton test particles of all N ens parallel ensembles are initialized 
according to Eq. (jSJ) with the same centroid (r^,p^), and the calculation is repeated for 
every event of the first stage. Thus, the antiproton test particle contributions to the mean 
mesonic fields reflect the presence of a real antiproton at the phase-space point (r^p^). 
In this new calculation, therefore, the compressional effects will manifest themselves in full 
strength without dilution. Further evolution of the p-nucleus system is calculated in a similar 
way as in 26j by using the GiBUU model without annihilation. However, in distinction to 
261 ]. we now take into account all collisional channels different from the annihilation one, 
in particular, NN — > NN and NN — > NN. This models dissipation leading to some 
small heating of the nuclear system and slowing down the antiproton during compression 
process. For brevity, sometimes we refer to the GiBUU calculations without annihilation as 
"coherent" calculations below. 

Instead of explicitly treating the annihilation on the second stage of calculations, we 
compute the survival probability of an antiproton in the course of compression as 

P surv (t) = exp j- J dt' r ann (0 1 . (6) 

Here 

Tann =< ^rel^anti > p (7) 

is the antiproton width with respect to the annihilation, p is the local nucleon density, v re \ 
is the relative velocity of an antiproton and a nucleon and cr ann is the p-annihilation cross 
section on a nucleon. Brackets in Eq. (J7J) denote averaging over the nucleon Fermi motion. 



The two-stage scheme described above is not fully equivalent to the true molecular dy- 
namics simulation. However, the most interesting phenomenon which we want to study, i.e. 
the dynamical compression of a nucleus by a slow antiproton, can be realistically simulated 
in this way. 

As we will see below (c.f. Figs. [Hand [7]), the width a r of the Gaussian in Eq. (J5J is a very 
important parameter, which can not be determined from our model. We, therefore, consider 
two most representative values: a r = 1 fm and o r = 0.14 fm. The first choice corresponds to 
a rather wide wave packet which presumably describes the static wave function of a strongly 



bound antiproton implanted in a nucleus 



17 



18 



26| . The second choice of a narrow wave 



packet is adjusted to describe the charge r.m.s. radius of a physical (anti)proton, r p = 0.9 
fm 34|. Indeed, in our model, the true source charge distribution of an antiproton is given 
by folding the coordinate space Gaussian ([5]) with the test particle Gaussian. Thus, we 



have a p = ^Ja^ + L 2 , where a p = r p /\/3 is the charge distribution width of a physical 
(anti)proton. 

The second-stage calculations can be significantly accelerated if one neglects the changes 
in a target nucleus caused by the antiproton cascade on the first stage. Then, the spherical 
symmetry of the 16 target nucleus can be utilized. In this case, the compressional evolution 
depends only on three variables: the absolute values of the antiproton initial radius-vector 
r and momentum p, and on the angle G = arccos(rp/rp) between r and p. (One needs six 
variables r, p in the case of arbitrary shape). Therefore, the second-stage calculations have 
been performed with the target nucleus for the set of the antiproton initial positions r and 
momenta p taken on the uniform 7 x 20 x 9 grid in the space (r, p, cos G), where r G [0.5; 3.5] 
fm, p G [0.05; 1.00] GeV/c, and cos G [—1; 1]. The results of the second-stage calculations, 
in particular, the antiproton survival probabilities at the time moments corresponding to the 
system entering to and exiting from the compressed state, have been stored. To determine 
the compression probability for a given first-stage event, resulting coordinates and momenta 
of the antibaryon at the annihilation point have been projected on the grid. 



III. DYNAMICS OF NUCLEAR COMPRESSION 

In this section, the nuclear response to the moving antiproton is considered disregarding 
p-annihilation. The latter is, however, implicitly taken into account by following the time 

8 



dependence of the p-survival probability. 

Figure [I] shows the time evolution of the nucleon and antiproton density distributions 
for the p l6 system for different p initializations. Only the cases are presented, when the 
initial antiproton momentum p is (anti)parallel to the initial position vector r, i.e. O = 
and O = 7r. 

Let us start by considering how compression depends on the initial antiproton coordinate 
z at fixed momentum p = 0.3 GeV/c. If the antiproton moves towards the nuclear centre, 
i.e. O = 7T, the compression of a nuclear system up to densities ~ 2p is reached within the 
time interval of about 10 fm/c (see panels (a)-(d) and (e)-(h) of Fig. [TJ. For the outgoing 
antiproton (O = 0), the compression is much smaller (panels (i)-(l) of Fig. [TJ, since the 
antiproton moves through the nuclear periphery. It is interesting, that at p — 0.3 GeV/c 
the antiproton does not leave the nucleus but only bounces off the nuclear boundary and 
finally gets captured. However, the capture takes place on the time scale of ~ 20 fm/c and, 
therefore, would hardly be observed due to a very low survival probability of the antiproton 
(see Fig. [2] below). 

For a higher momentum p = 0.8 GeV/c (panels (m)-(p) of Fig. [1]), the compression is 
practically absent, since the nuclear response is much slower then the time needed by the 
antiproton to cross the nucleus. We also see, that at p — 0.8 GeV/c the antiproton escapes 
from the nucleus, because its total in-medium energy Ep = \Jp 2 + mp 2 + Vp exceeds the 
vacuum mass by about 165 MeV. Here mp = + £(m* N — m^v) ~ 0.91rriN is the 
antiproton effective mass and Vp = — (308£) MeV ~ —68 MeV is the antiproton vector 
potential at p = p 

The compression process is quite sensitive to the choice of initial Gaussian width of the 
antiproton. One can see this from Fig. [T]by comparing the panels (e)-(h) and (q)-(t), where 
the calculations are shown for the same initial positions and momenta of p, but for the 
different widths a r . Due to a deeper nucleon potential, a smaller initial p-width makes the 
compression more pronounced and fast. Unless stated otherwise, the case of a r = 1 fm is 
discussed below. 

In Fig. [21 we present the time dependence of the maximum nucleon density p max and of 
the antiproton survival probability (|6]) for various antiproton initializations shown in Fig. CD 
As we have already seen in Fig. (TJ for the initializations with z < (i.e. for O = tt) and 
p = 0.3 GeV/c, nucleon densities up to 2p are reached. The antiproton survives with the 

9 
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FIG. 1. (Color online) Nucleon (thick solid lines) and antiproton (thin solid lines) densities as 
functions of the longitudinal coordinate z at different time moments for the p 16 system. The 
antiproton has been initialized on the axis passing through the nuclear centre, i.e. x = y = with 
momentum p along the positive z-direction. Different rows correspond to different p initializations 
characterized by the Gaussian width a r (fm), momentum p (GeV/c), and coordinate z (fm): (a)-(d) 
- (a r , p, z) =(1, 0.3, -2.5); (e)-(h) - (1, 0.3, -0.5); (i)-(l) - (1, 0.3, 2.5); (m)-(p) - (1, 0.8, -0.5); 
(q)-(t) — (0.14, 0.3, -0.5). The antibaryon annihilation is switched off in this calculation. 
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FIG. 2. (Color online) Maximum nucleon density (a) and antiproton survival probability (b) as 
functions of time for the p 16 system. Different curves correspond to different ^-initializations as 
explained in Fig. [TJ 
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probability P S urv ~ 10~ 2 until the time moment when the maximum density p max = 2po is 
achieved. 

The nuclear compression caused by an antiproton could only be observed, if the antiproton 
would annihilate in the compressed nuclear environment. This process can be detected by its 



specific final state characteristics. As shown in [18|, |26j], possible observable signals include 
the enhanced radial collective flow of nuclear fragments, hardening the energy spectra of 
emitted nucleons, and softening the meson invariant mass distributions. Moreover, the multi- 
nucleon annihilation (MNA) channels with the baryonic number B > 1 might be enhanced if 
the compressed zone is formed. A more exotic scenario, the deconfinement of an annihilation 
zone leading to the enhanced strangeness production has also being discussed in literature 
jisl . [35L 37]. Herein, we do not consider any specific signals caused by annihilation in the 
compressed nuclear state. We rather concentrate on the evaluation of the total p-annihilation 
probability at enhanced nucleon densities. For brevity, we refer to this possibility as to the 
annihilation in a compressed zone (ACZ) below. 

Let us define the compressed nuclear system as a system where the maximum nucleon 
density p max exceeds some critical value p c . If not stated otherwise, we choose p c = 2po in 
calculations. Such density values can be reached, e.g. in central heavy-ion collisions at beam 
energies of hundreds MeV/nucleon [381 ] . The probability for the antiproton to annihilate at 
Pmax > Pc is defined as 

-fann = ^surv(^l) — -Psurv^) , (8) 

where the time interval [tijt^] encloses the high-density phase of the time evolution, i.e. 
Pmax(ii) = Pmaxfe) = Pc with p max (t) > p c for t x < t < t 2 H For example, in the case 
(ov, p, z) =(1 fm, 0.3 GeV/c, -0.5 fm) we obtain t\ = 8.4 fm/c and t 2 = 11 fm/c (see Fig. [2]). 
Since the p survival probability drops exponentially with time, we have P surv (ti) 3> Psunih) 
and, therefore, actually P a „ n ~ -P S urv(^i)- 

Figure [3] shows the antiproton ACZ probability as a function of the p initial radial po- 
sition and momentum. As expected, the p-initializations with smaller momentum lead to 
larger P^ nn . The radial dependence of P£ nn at fixed initial momentum is somewhat more 
complicated. In the case of a larger width of the initial antiproton space distribution (<r r = 1 
fm), Pa nn has a weak maximum at r ~ 1 — 2 fm and decreases towards the nuclear centre 



When there arc more than one such intervals, the earliest one is chosen. 
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FIG. 3. (Color online) The contour plots of the antiproton ACZ probability P^ nn at p max > 2po 
(see Eq. ([8])) in the plane given by initial values of the antiproton radial position r and momentum 
p for the system p 16 0. The values of P£ nn are averaged over the cosine of the angle between the 
initial radius vector and momentum of the antiproton. Panel (a) ((b)) corresponds to the initial 
antiproton width a T = 1 (0.14) fm. 
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slightly. This can be traced back to Fig. [TJ where we see, that the p initializations at different 
positions result in practically the same compressional effect provided that the antiproton 
moves to the nuclear centre (c.f. panels (a)-(d) and (e)-(h)). For a narrower initial antipro- 
ton space distribution (a r = 0.14 fm), the maximum of the ACZ probability is located at 
the nuclear centre, since compression is much faster in this case, and, thus, is more sensitive 
to the local nucleon density. 



IV. p-NUCLEUS COLLISIONS 

As it was demonstrated in the previous section (c.f. Fig. [3]), the ACZ probability depends 
on the position and momentum of the antiproton at the beginning of compression process. 
Therefore, before discussing the results of a full two-stage calculation, it is instructive to 
study the distributions of antibaryon annihilation points in the coordinate and momentum 
space. These distributions are determined at the first stage of calculations. 

Figure H] shows the radial distributions of the antibaryon annihilation points for the 
p 16 reaction at several beam momenta. For inclusive events (a), the maxima of these 
distributions are located at the peripheral region, where the density is about 30% of the 
central density, independent of the beam momentum. This is a pure geometrical effect caused 
by mixing of events with all possible impact parameters. However, for central collisions (b), 
the maxima are shifted closer to the nuclear centre. The shift becomes larger at higher 
beam momenta. This is expected, since with increasing pi a b the p-nucleon annihilation cross 



section drops [39| leading to the larger fraction of deeply-located annihilations. 

Figure [5] demonstrates the momentum distributions of antibaryons at the annihilation 
points. There is only a little difference between the shapes of the distributions for the 
inclusive (a) and central (b) events at the same beam momentum. However, the total 
annihilation probability is increased by a factor of 3-10 for the central collisions, which is 
also seen in Fig. HJ The distributions have a sharp peak at the beam momentum and a long 
tail towards small momenta. In the case of the smallest beam momentum p^h = 0.3 GeV/c, 
the peak is broader and shifted to the higher momenta p > pi^ for the central collisions. 
This is caused by the antiproton elastic collisions with the Fermi sea nucleons and by the 
antibaryon acceleration in a strongly attractive mean field potential. 

The acceleration is better visible in Fig. [6] which shows the correlation between the radial 
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FIG. 4. (Color online) Radial distributions of annihilation points for p 16 collisions at differ- 
ent beam momenta (see key notations) normalized to the total annihilation probability. For the 
reference, thin solid lines show the nucleon density profile of the 16 O nucleus. Panels (a) and 
(b) represent the inclusive spectra (b < 8 fm) and the spectra for central collisions (b < 1 fm), 
respectively. Note different scales of vertical axes in panels (a) and (b). 
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FIG. 5. (Color online) Momentum distributions of annihilation points for p 16 collisions at different 
beam momenta (see key notations). Panels (a) and (b) represent the inclusive spectra and the 
spectra for the central collisions, respectively. 
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FIG. 6. (Color online) Average momentum of annihilating baryon as a function of the radial 
position for the inclusive set of events p 16 at different beam momenta. The error bars represent 
the dispersion of momentum distribution at a given r. 

position and momentum of an antibaryon at the annihilation point. The centrality depen- 
dence is quite weak in this case, thus we have presented the results for the inclusive event set 
only. For pi a b = 0.3 GeV/c, the average momentum of annihilating antibaryon increases up 
to 0.5 GeV/c at the nuclear centre. For larger beam momenta, the mean field acceleration 
is hindered by the collisional damping of an initial p momentum. 

We will now discuss the results of the full two-stage calculations (see Sect. [TTJ) . The total 
annihilation cross section on a nucleus cx ann and the ACZ cross section a corapr are determined 
as follows: 

-1 iV an n(&) 



a, 



£ 2VT&A6-— - J2 P a c nn (r,, Pl ) . (10) 



compr " ^ N (b) 



Here, N ev (b) and A r ann (6) are, respectively, the total number of events and the number of 
annihilation events calculated within standard GiBUU (the first stage) for a given impact 
parameter b. The quantity P a ^ n (rj, p 4 ) (see Eq. (jS])), which depends on the annihilation point 
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position rj and momentum pj in a given annihilation event i, is the annihilation probability 
at p max > p c computed within a coherent GiBUU run (the second stage). The cutoff value 
of the impact parameter 6 max has been chosen to be 8 fm for an inclusive event set and 1 fm 
for central events. 

Figure [7| shows the annihilation cross section o" ann (a) and the relative fraction of ACZ 
Ccompr/cann (b) as functions of the beam momentum. While a ann drops with increasing pi a b 
due to the momentum dependence of the pN annihilation cross section, the ratio o" compr /cr ann 
reveals an interesting nonmonotonic behaviour. First, it drops with increasing beam mo- 
mentum up to pi a b — 1 GeV/c and then starts to increase saturating at pi a b ~ 3 GeV/c. The 
growth of this ratio at pi a b > 1 GeV/c is caused by opening the inelastic production chan- 
nels, NN ->■ NNtc with the threshold beam momentum p thr = 0.787 GeV/c, NN ->■ NNtttc 
with pthr = 1-210 GeV/c etc. The inelastic production leads to the additional deceleration 
of an antibaryon and, therefore, increases the probability of the nuclear compression [l^ 
(see also Fig. [S]). Selecting the central events increases the ratio Oc^^/a^ by about a 
factor of three, which is caused by a larger relative fraction of annihilations at small radii 
(c.f. Fig. [lb). 

The important result of the previous section is that only a slow antiproton can induce nu- 
clear compression. In practice, we have used the ensemble of annihilation points to initialize 
the coherent GiBUU runs assuming that antiprotons become slow not far away from their 
annihilation points. To check this assumption, we have performed additional calculations 
with other transition criteria from collisional deceleration stage to the coherent compression 
dynamics. In the first calculation, we have generated the ensemble of points where the 
momenta p and coordinates r of antibaryons satisfy the conditions |p| < p c and |r| < r c 
simultaneously, i.e. when the antibaryons become slow enough and close enough to the nu- 
clear centre. Here, p c and r c are parameters to be chosen. As follows from Fig. [31 the choice 
p c ~ 0.3 — 0.5 GeV/c, r c ~ 3 fm provides almost the full coverage of the (r, p)-region where 
a significant (p > 2p ) compression is expected. Another criterion selects the antibaryon 
momentum and position at the first time instant, when the antibaryon becomes bound, i.e. 
its energy falls below its bare mass. Figure [8] shows the ACZ probability calculated by using 
the different transition criteria. All results are quite similar, except for the calculation with 
p c = 0.5 GeV/c which becomes unphysical at pi ab < p c . 

To give more insight into the p-induced nuclear compression, in Figs. [9] and [10] we show 
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FIG. 7. (Color online) Antiproton annihilation cross section on 16 (a) and the probability of 
annihilation at /0 max > 2po (b) vs a beam momentum. Different lines correspond to different 
centralities: b < 8 fm — inclusive set of events, b < 1 fm — central events. The cases of a narrow 
(cr r = 0.14 fm) and a wide (ay = 1 fm) initial antiproton space distribution are depicted in panel 
(b) by lines with open and filled squares, respectively. 
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FIG. 8. (Color online) Same as in Fig. [7J (b) for the inclusive event set with oy = 1 fm, but 
for different criteria of transition to the coherent compression dynamics (see text). The lines with 
filled and open circles represent, respectively, the results obtained with a criterion using the critical 
momentum p c = 0.3 GeV/c and 0.5 GeV/c, respectively. The line with filled triangles corresponds 
to the criterion requiring that the antibaryon becomes bound. The line with filled squares shows 
the calculation with the default criterion using the annihilation points, same as in Fig. [7] (b). 
Additionally, the line with open pentagons shows the results obtained by switching-off the inelastic 
channels of the NN scattering. In this case, the ACZ probability quickly drops with increasing 
beam momentum and becomes less than 10~ 8 at pi a b > 1 GeV/c. 

the distributions of annihilation events on the longitudinal coordinate z and the longitudi- 
nal momentum p z of the antibaryon for central (a) and peripheral (c) p w O collisions. The 
relative probability of annihilation in the compressed zone as a function of z and p z is shown 
in panels (b) and (d) for central and peripheral collisions, respectively. Independent of the 
beam momentum, the maximum ACZ probability is reached if the antibaryon is stopped 
in the central region. However, the longitudinal coordinates for events most favourable for 
compression are rather uncertain, as expected already from Fig. [TJ On the other hand, we 
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FIG. 9. (Color online) Panels (a) and (c) show the contour lines of the antiproton annihilation cross 
section d<7 ann / (dzdp z ) (mb/(fm • GeV/c)) in the plane of longitudinal coordinate z and longitudinal 
momentum p z of the annihilation points for central (a) and peripheral (c) collisions. Panels (b) 
and (d) show the contour lines of relative fraction of annihilations at high density (/9 m ax 
<7comp/c"ann, in the same plane for central (b) and peripheral (d) collisions. The colliding system is 
p w O at 0.3 GeV/c. 

definitely observe a rather strong impact parameter dependence with the clear preference of 
central collisions for selecting the ACZ events. At large beam momenta (Fig. [TO]) . the com- 
pression can only be reached at the extreme tail of the antibaryon longitudinal momentum 
distribution, and the total probability of ACZ is small. As one can see from the right panels 
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FIG. 10. (Color online) Same as Fig. [H but for pi a b = 3 GeV/c. Notice different scales of vertical 
axes in the left and right panels. 

in Figs. M and [101 a significant compression (p max > 2p ) can be produced by antibaryons 
whose longitudinal momenta are less than 200 MeV/c. However, the maximum relative 
probability of ACZ in the (z,p z )-pla.ne is practically independent on the beam momentum. 
This is expected, since in our model the probability of ACZ depends only on the position 
and momentum of the antibaryon prior the annihilation. 

Since nuclear compression is most probable for stopped annihilations, one needs a trigger 
to select the events with slow antiprotons. We will discuss two possible triggers here. 

The first trigger requires the detection of a fast proton in forward direction 



2JJ. The 
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idea behind is that the incoming antiproton can be decelerated and captured in a nucleus 
by experiencing a hard collision with a single nucleon. Figure [IT] (a) shows the cross section 
<7 Pmin of an antiproton annihilation on 16 O accompanied by the emission of a proton with 
momentum exceeding some value p m j n as a function of p m i n . For simplicity, we did not apply 
any angular cuts for the proton momentum. At large beam momenta, 3 and 10 GeV/c, the 
cross section cr Pmin sharply drops with p min near p min ~ pi^. In the panel (b) of Fig. [II], 
we show the relative fraction cr compr /o" Pm . n of ACZ as a function of the minimum proton 
momentum p m i n . For p^ = 3 and 10 GeV/c, the quantity a compr / 'erp min grows by almost a 
factor of thirty while p m i n increases from zero to pi^. 

Emission of a fast proton with momentum close to the p-beam momentum can be caused 
by the following mechanisms: (i) elastic scattering p + p — > p s i ow + Pf ast , (ii) inelastic pro- 
duction processes of the type p + p — > p s \ ow + pf ast + mesons, and (iii) collisions with high- 
momentum annihilational pions, n+p — > Pf as t + X. We have checked, that inelastic reactions 
(ii) give the largest contribution to the production of the fast proton at pi a b = 3 and 10 
GeV/c. This makes the fast proton trigger rather efficient at high beam momenta. The 
contribution from process (iii) is relatively small, while elastic scattering (i) practically does 
not contribute to the yield of fast protons at pi a b = 3 and 10 GeV/c. On the other hand, for 
small beam momenta, 0.3 and 1 GeV/c, the pionic mechanism (iii) contributes dominantly 
to the fast proton yield with only a small admixture of elastic scattering (i). Therefore, the 
trigger based on a fast proton in final state is ineffective at small beam momenta. 



The second possible trigger is based on the energy deposition [40|, |41| 



E dep = T p + 2m N -Y,E£ es , (11) 



where T p is the antiproton beam energy, E$ es is the energy of i-th. outgoing meson, and the 
sum runs over all produced mesons in a given event. Neglecting nucleon and antibaryon 
binding energy, antiproton elastic and inelastic scattering before annihilation and final state 
interactions of produced mesons, one has -Edep = 0. In the case of low energy antiproton- 
nucleus collisions, annihilations with a larger energy deposition occur deeper in the nucleus, 
as was found in [40] . The explanation was that the annihilation mesons loose their energy or 
get absorbed more effectively if the annihilation takes place deeply inside the nucleus. For 
high-energy p-nucleus interactions, the incoming antiproton can rescatter before annihilation 
transferring a part of energy to the nucleons. This also leads to larger values of E dep , since 
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FIG. 11. (Color online) (a) The cross section cr Pinin of p- annihilation on 16 O accompanied by the 
emission of a proton with momentum larger than p m i n as a function of p m i n . (b) The relative 
probability cr coiapT /a Pm . m of the annihilation in compressed zone (/9 max > 2po) vs the minimum 
momentum p m i n of the emitted proton. Different curves refer to different p-beam momenta. 
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the produced mesons will have a smaller total energy in this case. Both types of events, 
deep and/or slow antibaryon annihilations, should be of the ACZ-type with an increased 
probability. Results for i?d ep -trigger are shown in Fig. [T2J As one can see, triggering on a 
large energy deposition, E^ep — Tp + 2rriN, increases the fraction of ACZ events by about a 
factor of thirty with respect to E^ cp ~ for the beam momentum of 10 GeV/c. 



V. IN MEDIUM MODIFICATIONS OF ANTIPROTON ANNIHILATION 



So far we have used the vacuum pN annihilation cross section and the fixed antibaryon 
mean field parameters. As the survival probability of an antiproton ([6]) is determined by 
its annihilation width (J7J), it is important to consider possible in-medium modifications of 
the latter. On the other hand, the speed and the amplitude of the nuclear compression 
depends on the value of the reduction factor £ of the antibaryon-meson coupling constants 



261 ] . Thus, the probability of ACZ is the result of the competition between compression and 



annihilation. In this section, we discuss possible modifications of the antiproton annihilation 
width in nuclear medium. The sensitivity of our results to the in-medium modifications is 
demonsrated in Fig. [15] below. 

As discussed by many authors (see e.g. Refs. 



1* 



23 



25 



35 



42 



46| ) , the annihilation rate 



of antiprotons in a dense nuclear medium may significantly differ from simple calculations 
using the vacuum pN annihilation cross section o~ ann (see Eq. (JTJ)). There are several effects 
which become important at sufficiently high nucleon densities. First, the dispersion relations 
of nucleons and antinucleons are modified due to interactions with mean mesonic fields. In 
particular, the effective masses m* N and m*^ are reduced compared to the vacuum value. 
Generally, this leads to the reduction of the imaginary parts of the nucleon and antinucleon 
self-energies in nuclear medium. The influence of the baryon and antibaryon in-medium 

. Other 



dispersion relations on the antibaryon annihilation rate has been studied in Ref. [18] 
examples of the influence of the baryonic effective masses on hadronic processes are the in- 
medium reduced cross sections of the NN elastic scattering 47], |48j and of the A-resonance 
production NN — > NA 49|, [50[] . As an illustrative example of the in- medium reduction 
caused by effective masses, we consider the two-pion annihilation channel later-on in this 
section. 

Another important in-medium effect is the appearance of the MNA channels. The famous 
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Pontecorvo reaction 5l| pd — » vr p is an example of the MNA processes. It is commonly 



believed that MNA is responsib 
p annihilation on nuclei 



24 



52 



e for the emission of high-energy protons from low-energetic 



531 ] . Moreover, the triggering on high- momentum protons 



is already applied in experimental techniques to distinguish MNA from the single-nucleon 
annihilation (SNA) followed by the final state interaction (rescattering and absorption) of 
produced mesons [53 ]. 

Up to now the attempts to estimate the MNA contribution performed by different the- 
oretical and experimental groups did not lead to definite conclusions. The experimental 
determinations of the MNA probability give the values of about 10-30% for the p annihila- 
tions at rest 52|, 53 1. One has to admit that these values agree with predictions of Hernandez 



and Oset (HO) 23j-|25j. It is important for this agreement, however, that the annihilations 
of stopped antiprotons take place at the nuclear fringe, p ~ O.lpo 52j. On the other hand, 
HO argue in Ref. 24J that the MNA channels are required to describe the high-momentum 
tails of the proton momentum spectra from p annihilation at pi a b = 608 MeV/c on 12 C 54]. 



and the GiBUU calculations 



22| have 



However, the intranuclear cascade calculations 

demonstrated that the agreement with the data [54j can be reached without any unusual 
mechanisms. 



Using a diagram language, HO have considered pN annihi 



ation vertices including vir- 



tual pions which may decay into particle-hole excitations [23j-|25j. These diagrams can be 
interpreted as MNA channels, which, according to the HO calculations, have extremely high 
probability at p ~ po, one order of magnitude higher than the ordinary SNA channels. This 
result is in a clear contradiction with the theoretical estimates by Cugnon and Vandermeulen 
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43] and by Mishustin et al. 18], although these estimates are based on relatively simple 



statistical considerations. In our opinion, the HO calculations have a problem with conver- 
gence of the series in powers of p at p > p Q (Eq. (4.3) in Ref. [24J). Since it is very difficult 
to test the MNA probability at p ~ po experimentally, different theoretical predictions are 
still possible here. 

In Ref. 18], the relative importance of MNA channels was evaluated by calculating the 
probability to find more than one nucleon in the annihilation volume V ann . This calcula- 
tion was done for a spherical volume with the radius R ann — 0.8 fm assuming the Poisson 
distribution in the number of nucleons n, P(n) = n n exp(— n)/n\, where n = pV ann is the 
average number of nucleons in this volume. In the case of enhanced density, p ~ 2p , one 
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has n ~ 0.6. This leads to the probability of MNA channels with n > 1 on the level of 40% 



of the SNA 
in Refs. 23J 



n = 1), which is about one order of magnitude smaller than the value predicted 



25|. 



In the literature, one can also find other arguments against a large enhancement of pN 



annihilation cross section in nuclear medium. For instance, as argued in Refs. |44|, |45j, 
the presence of additional nucleon(s) in the annihilation volume may lead to the "screen- 
ing" of in-medium annihilation as compared to the usual SNA mechanism. By introducing 
the screening effect these authors achieve a better agreement with experimental data on p 
production in pA and AA collisions at AGS energies. 

In order to illustrate the influence of the in-medium effective masses on the antiproton 
annihilation, let us consider a relatively simple case of the two-pion annihilation pp — > 7t~tt + . 



56 



57]. In 



In vacuum, this process can be described by the one-nucleon exchange model 
the Born approximation, the matrix element can be written as follows: 

M = -2F 2 (t)v(m sp ,pp)[A- fcB(t)]u(m sp ,p p ) , (12) 

where u and v are, respectively, the proton and antiproton bispinors (uu — 1, vv — —1), 
which depend on the spin magnetic quantum numbers m sp , m sp and on the four-momenta 
p p , p.p, and k is the four-momentum of tt + . The scalar parameters A and B are defined as 

A=^2m N , 5(t) = 4(l+7^) . (13) 

where / = 1.008 is the pion-nucleon coupling constant and t = (p p — k) 2 . The factor of -2 
in 



56j: 



121) is obtained from an isospin algebra. The off-shell nucleon form factor is chosen as in 



where A is a cut-off parameter. In the center-of-mass (cm.) frame, the differential cross 
section of the process pp — > 7i~n + is 

d<Jpp^ 1T - 1T + _ (2m jy) i^|2 ^ 7r7r (15) 
dVL c , m , 647r 2 s q Pp 

where s = (p p + p p ) 2 , qp p = q(\fs,rriN) and q n7r = q(y/s,m n ) are the cm. momenta of 
the incoming and outgoing particles, respectively, with q(y/s,m) = (s/4 — m 2 ) 1//2 , and 



\M\ 2 — j \M\ 2 . Solid line in Fig. Unreports the total vacuum pp — > n n + cross section 
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FIG. 13. The cross section of the process pp — > tt~tt + as a function of the total in-medium cm. 
kinetic energy y/s* — 2m* N . The results are shown for the different choices of a nucleon effective 
mass m* N . The calculation with m* N = mj\T (thick solid line) corresponds to the vacuum cross 
section, which has to be compared with experimental data. The calculations with m* N = O.Gtotv 
(thin solid line — full result, dashed line — vacuum matrix element) and m* N = 0.3m^ (dotted line 
— full result, dash-dotted line — vacuum matrix element) represent the in-medium cross sections. 
Experimental data are from Refs. 58l462[|. 
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calculated in the Born approximation. This is a quite rough approximation in the case of pp 
incoming channel. We stress, however, that our main purpose here is just to demonstrate 
the influence of the in-medium effects and not to perform the state-of-art calculations for 
the two-pion annihilation in vacuum. To fit the data for slow antiprotons (p\ a b < 1 GeV/c), 
we have chosen a rather small value of the cut-off parameter A = 1.0 GeV, since the initial 



state interactions are neglected (see discussion in Refs. [56l. l57j|). 

Assuming for simplicity the G-parity transformed proton scalar and vector potentials 
acting on the antiproton, the baryonic mean fields can be now taken into account by replacing 
m N m* N , p p -»■ p*, PjL-tP* p , s* = (p*+p* p ) 2 , and t t* = (p* - k) 2 in Eqs. , QS) 



and (IT5|) (c.f. Refs. 50|, l63j). Note, that we always keep the vacuum nucleon mass in the 
numerator A 2 — m 2 N of the form factor (|14[) . since the above replacements are motivated by 
the baryon in-medium self-energies which should not change the form factor fixed in vacuum. 
Then, the total in-medium pp — » 7r~7r + cross section reads: 



^med , .nzx _ (2rn* N ) 2 q(^,m w ) 



1 

J dcos 9 c . m . |M mcd | 2 (VF, cos 9 c . m .) . (16) 



327T s* q(\fs* , m* N ) J 

Here, the quantity M med is the in-medium matrix element, while c . m . is the angle between 
the proton and 7r + three-momenta in the cm. frame. The results of calculation using Eq. 
( TIE]) at m* N = 0.6mjv (p = po) and m* N = 0.3rriN (p — 2p ) are shown in Fig. [13] by the thin 
solid and dotted lines, respectively. As one can see, the in-medium pp — > cross section 

is strongly reduced, largely due to the (2m* N ) 2 factor in Eq. (TIB"]) , which comes from the 
in-medium Dirac bispinor normalization. For orientation, we present in Fig. [TH the baryon 
density dependence of the nucleon effective mass m* N and of the nucleon scalar density ps 



see Eqs. (8), (9) and (15) in Ref. 26[)- The effective mass m* N drops with increasing baryon 



density which is an important effect influencing in-medium cross sections (see also Refs. 



47H50J.) 

The two-pion annihilation channels represent, however, less than 1% of the total annihi- 
lation cross section, while the direct calculation of multi-meson annihilation matrix elements 
is impossible. We will assume that the matrix elements are not modified in nuclear medium 
and take into account only the in-medium bispinor normalization, flux and phase space vol- 
ume. This intuitive assumption has some support from the earlier studies of NN — > NN 
and A^A^ — > NA cross sections in nuclear matter 



47 



50} . In this way, one can write the 



following formula for the in- medium cross section of the AW — > Mi, M 2 , M n annihilation 
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FIG. 14. (Color online) Nucleon effective mass and scalar density vs the baryon density in infinite 
nuclear matter in the case of NL3 interaction [30] applied in the present work. 



channel with n outgoing mesons (c.f. 26|, |28|, |64j): 



mcd 



fm) 



N \ 1 NN 



"A \ - •.!/, . U; M„( V S*) — 0NN-+Mi,M 2 ,...,M n (V s corr) ( " , ? 

$ n (A/g*; mi,m 2 , -,rn w ) 
corn mi,m 2 , ...,m n ) 

Here m 1 ,rri2, ■■■,m n are the vacuum masses of outgoing mesons, i/^Wr" = — 2(m^ — 
m^v) is the so-called corrected invariant energy of the colliding particles, the analogue of 
the invariant energy in vacuum. The quantities I^ N = q ( y/s corr ,m N ) ^s corr and I^ N = 
q{\fs* ) m* N )\fs* are the vacuum and in-medium flux factors, respectively. The n-body phase 
space volume is defined as 

d 3 k x r d 3 k 2 f d 3 k r . 



$ n (v^;m 1 ,m 2 ,...,m^ = /-^i-/ 



(2tt) 3 2u 1 J (2n) 3 2u 2 
x S {4) (V-k 1 -k 2 -...-k n ) , 



(27r) 3 2o; r , 



(18) 



where V 2 = s and uf — 



mf U 



1,2, ...,n). The quantity ^^ Mi ,m 2 ,...,m„(Vw) 
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is the vacuum cross section. Application of Eq. (JT7J) to the process pp — > it 7r + leads to 



the formula (USD with replacement \M med \ 2 (y/s*, cos 9 c . m .) \M\ 2 (y/s^, cos c . m .). The 
results obtained using Eq. (fT7|) for = 0.6m N and = 0.3m7v are shown in Fig. [TBI by 
the dashed and dash-dotted lines, respectively. The conclusion is that using vacuum matrix 
element produces somewhat less pronounced in-medium reduction of the cross section. 

In order to simulate the mean field effects on the total annihilation cross section, we 
represent it as a sum of partial annihilation cross sections with various outgoing mesonic 



channels. In practice, this is done by using the statistical annihilation model of Refs. [651.166] . 
Then we apply Eq. (fT7|) to every annihilation channel with up to n = 6 outgoing mesons. 
The mean field effects on the annihilation channels with more than six outgoing mesons are 
neglected. 

In the following GiBUU calculations of the present section, we keep the first stage (see 
Sect. [TI]) unchanged: it is always computed with the reduction factor £ = 0.22 and vacuum 
pN annihilation cross sections. This is reasonable, since the compressional response of a 
nuclear system on a fast moving antiproton is weak and can be neglected. On the other 
hand, the in-medium corrections to the pN annihilation cross section should also weaken for 
the fast antiproton. Thus, we vary the model parameters for the second stage compressional 
dynamics only. 

Fig. [15] shows the probability of p annihilation at p max > p as a function of p for the 
inclusive set of p 16 events at the beam momentum of 3 GeV/c. First, we study the influence 
of the reduction factor £ of antiproton-meson coupling constants on the ACZ probability. To 
this aim, we have performed additional calculations by choosing £ = 0.15 (Re(K pt ) = —105 
MeV) and £ = 1 (Re(V opt ) = —677 MeV), where the values of the Shrodinger equivalent 
potential at = (see Ref. [22J for details) in the centre of the 16 nucleus are given in 
brackets. We recau that o Ur default value £ = 022 (Re(K P .) = -153 MeV) is motl vated by 
the best agreement of the GiBUU calculations [22[ with the measured antiproton absorption 
cross sections on nuclei. The value £ = 0.15 is in line with the results of p- atomic X- 
ray transitions and radiochemical data analysis [19j, while the extreme choice of £ = 1 
corresponds to the G-parity transformed nucleon mean fields. As expected, larger (smaller) 
values of £ give rise to larger (smaller) ACZ probability at a given p. Results are quite 
sensitive to the antiproton-meson coupling constants. For instance, at p = 2po, the ACZ 
probability is equal to zero in the case of £ = 0.15, i.e. the maximum nucleon density does 
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FIG. 15. (Color online) Probability of p annihilation at p max > p vs p for p 16 collisions at pi a b = 3 
GeV/c. The calculations with various values of the reduction factor £ and vacuum pN annihilation 
cross section are shown by lines with points denoted by the values of £ only. The calculation with 
the in-medium enhanced pN annihilation probability by the HO formula (Eq. (5) in Ref. 25J with 



parameters taken at the p kinetic energy of 50 MeV) and £ = 1 are shown by the dashed line with 
filled boxes. The result taking into account both the HO formula and the mean field reduction of 
the pN annihilation cross section according to Eq. (|17p is shown by the dashed line with upside- 
down filled triangles. The calculations for p < po are not shown, since they are influenced by the 
finite size of the (r,p, cos O) grid in the space of the antiproton radial position r, momentum p and 
cos = rp/rp (see Sect. HI]) . 

not reach the value of 2p at all with this value of £. At the same time, for £ = 1, the ACZ 
probability is ~ 10~ 2 at p — 2p which is three orders of magnitude larger than for £ = 0.22. 

Finally, we discuss the sensitivity of our results to the choice of the in-medium annihilation 
cross section related to the imaginary part of the antiproton optical potential 

Im(V opt ) = - l - < v rel aZ d > P ■ (19) 
Here, cx£^ d is the total in-medium pN cross section, which includes both (in-medium) an- 
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nihilative and Pauli-blocked nonannihilative contributions. Note, that in distinction to Eq. 
([7]) for the annihilation width, Eq. ffl9|) contains the total pN cross section. Applying Eq. 
(fn?|) for the case of the HO formula for the annihilation probability per unit length (see Eq. 



(5) in Ref. 25]) leads to an extremely deep imaginary part of the antiproton optical poten- 
tial, Im(V^pt) — —(1200 — 1700) MeV at the centre of the 16 nucleus, with the uncertainty 
due to the antiproton mean field. This is more than one order of magnitude larger than the 
value lm(V opt ) — —107 MeV in our default choice of model parameters j^j. 

At this point, certainly, one wishes to get some phenomenological estimates of the 
antiproton-nucleus potential depths. Unfortunately, it is very hard to get stringent phe- 
nomenological constraints on the optical potential of a hadron at the nuclear centre [671] . As 
it is known for a long time from pionic atoms and low-energy pion nucleus scattering, the 
different density shapes of the potential give the same result for the atomic and scattering 



data, while they strongly differ at p = Pn [68|, |69|. In a similar way, the p-nucleus scattering 
and absorption cross sections (see [70h72|| and refs. therein) and the p-atomic data analysis 
19l . 1671 . 1 731 ] result in quite uncertain real and imaginary parts of the p optical potential at the 
nuclear centre lm(V opt ) = -(70 - 150) MeV and Re(V opt ) = -(0 - 100) MeV. We stress that 
the actual uncertainty in the potential depths may be much more due to the extrapolation 
from far periphery of a nucleus using some assumed relation between the nuclear density 
and potential. Nevertheless, the known phenomenological values are in a fair agreement 
with our default model inputs. 

As one can see from Fig. [15], using the HO formula leads to eight orders of magnitude 
smaller ACZ probability at p = 2p as compared to the calculation with the vacuum anni- 
hilation cross section (see the lines with filled boxes and with filled diamonds). This is not 



surprising, since the p-dependent terms in Eq. (5) of Ref. 25( strongly enhance the annihi- 
lation rate at high densities. However, as discussed above, the mean field and phase-space 
effects should reduce the annihilation rate. 

Now we implement both effects simultaneously by introducing corresponding multiplica- 
tive factors to the vacuum pN annihilation cross section. The resulting ACZ probability 
increases by eight orders of magnitude with respect to the one given by the HO effect alone 
(see the lines with upside-down filled triangles and with filled boxes in Fig. IT5"j) . i.e. practi- 
cally brings it back to the original calculation with the vacuum annihilation cross section. 
Certainly, this is only a rough estimate of the in-medium effects in the annihilation cross 
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section. In our opinion, the full calculation in the spirit of Refs. |23j-|25|, but taking into 
account, in- addition, the baryonic mean fields is needed in order to obtain the realistic values 
for the antiproton width at high nucleon densities. 

In spite of large uncertainties in the in-medium properties of antiproton, we think that 
our standard choice of the model parameters, i.e. the vacuum pN annihilation cross section 
and the reduction factor £ = 0.22, is quite reasonable for the present study of compressional 



effects. As it has been shown within the GiBUU model in [22J, this set of parameters accounts 
for the p absorption data on nuclei at pi a b < 1 GeV/c and pion and proton production data 
from p annihilation on nuclei at 608 MeV/c. 



VI. SUMMARY AND OUTLOOK 



We have generalized our 



arevious analysis of the nuclear compression dynamics induced 



by an antiproton at rest 26j to the case of a moving antiproton. The p-nucleus collisions 



}eam momenta of 0.3 — 10 GeV/c have been simulated within the transport GiBUU 



at the 

model 27j with relativistic mean fields. In our two-stage calculational scheme, we apply, first, 
the standard parallel ensemble mode of GiBUU to determine the antibaryon coordinates and 
momenta at the annihilation point. We have studied in-detail the coordinate and momentum 
distributions of annihilation points at different beam momenta. This calculation is performed 
in order to evaluate the probability that the antibaryon has been slowed down and reached 
the nuclear interior before annihilation. Those rare events which satisfy these conditions are 
used as the input for a more detailed calculation. Namely, we perform the coherent GiBUU 

n 

runs [26J initializing the antiproton at the given momentum and position inside the nucleus 
and following the evolution of the p-nucleus system. In the coherent mode, the antibaryon- 
nucleon annihilation channels are switched off, but, instead, the survival probability of 
the antiproton is determined as a function of time. This allows to trace the compression 
process of the p-nucleus system in time and determine the probability of p-annihilation in 
the compressed nuclear configuration with the maximum nuclear density p max > 2p - 

The results of our study are quite sensitive to the actual strengths of the real and imagi- 
nary parts of the antiproton optical potential. E.g., by choosing Re(V^ pt (po)) — —100 MeV 
instead of our default Re(V op t(po)) — —150 MeV reduces the ACZ probability by two or- 
ders of magnitude. The -100 MeV value of the real part of antiproton optical potential is 
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consistent with the X-ray data from antiprotonic atoms and radiochemical data [l9j, while 
-150 MeV value is favoured by GiBUU calculations of the antiproton absorption cross sec- 
ions on nuclei and of the pion and proton momentum spectra from p annihilation on nuclei 



22( | . Another important source of uncertainty is given by the value of lm(V opt {po)) ■ In 



our standard calculations, we adopt Im(V^ pt (po)) — —100 MeV, which follows from a sim 



pie ^approximation and is consistent with Re(T^, pt (po)) — —150 MeV as shown in [22]. 



On the other hand, according to the model of Hernandez and Oset [23N25|. the antipro- 
ton annihilation rate is increased by about one order of magnitude with respect to the 
simple tp-approximation, even at the normal nuclear matter density. This will result in 
Im(V^ pt (po)) — —1500 MeV, which is well beyond the phenomenological expectations. If 
this were, indeed, the case, the ACZ probability would be 5-8 orders of magnitude smaller 
than in our standard calculations. 

With all these reservations in mind, we now summarize the results of our standard calcu- 
lations which use the phenomenological input parameters for the antiproton-nucleus interac- 
tion. In general, antiproton initializations in a central nuclear region with momenta of less or 
about the nucleon Fermi momentum lead to the maximum probability of annihilation in the 
compressed zone of the order of 10~ 3 — 10 _1 . The uncertainty is caused by unknown spatial 
spread of the antiproton distribution function. When combined with the actual antibaryon 
positions and momenta at the annihilation points determined from the first stage GiBUU 
simulation, this results in the ACZ probability ~ 10~ 5 — 10~ 3 for the beam momenta of 3 — 10 
GeV/c. We have found that, within this beam momentum range, the excitation function 
of the ACZ probability is very flat (c.f. Figs. [5]). Therefore, the range pi a b = 3—10 
GeV/c is quite well suited for the study of compressed nuclear systems. The beam momenta 
of about 1 GeV/c are clearly disfavoured, since the antiproton is not decelerated enough 
due to the smallness of the NN inelastic production cross section. At < 1 GeV/c, the 
ACZ probability grows up with decreasing beam momentum. However, additional triggers 



demanding a fast proton 2l| or large energy deposition 40|, |4l| are not very efficient for ACZ 
selection at small beam momenta. On the other hand, we have found, that these triggers 
increase the ACZ probability by more than one order of magnitude in the beam momentum 
range of 3-10 GeV/c. Such antiproton beams will be available at FAIR which would be the 
ideal place to search for the nuclear compression effects induced by antibaryons. By taking 



the expected luminosity L = 2 ■ 10 32 cm 2 s 1 for the PANDA experiment at FAIR 
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2l|, the 



ACZ rate can be estimated as Y = cr compr L ~ 10 2 — 10 3 s _1 , where <7 compr ~ 10~ 3 — 10~ 2 
mb is the ACZ cross section above 1 GeV/c (see Fig. [7]). Here, we would like to stress once 
again that, due to the presently not well known antiproton optical potential at p > po an d 
due to uncertain spatial spread of the antiproton distribution function, the above estimate 
of the ACZ rate has a rather large uncertainty. 

We have also shown that the selection of small impact parameter events increases the 
ACZ probability by a factor of 2-3. This selection could be reached, e.g. by triggering on 
the events with a small azimuthal asymmetry of secondary particles. 

Some signals associated with the ACZ events have already been discussed in Refs 



18 



26). 



But, unfortunately, no unique signal suggested so far can alone be sufficient to identify 
the nuclear compression unambiguously. Therefore, we believe that the combination of 
different signals, e.g. emission of a fast proton plus large collective flow energy of the nuclear 
fragments, would be a more promising strategy to search for the ACZ events. Certainly, 
further theoretical studies are needed in order to find the experimentally realizable ways to 
observe nuclear compression in p-nucleus collisions, in particular at FAIR energies. 
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